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Abstract. Using subspace methods, we study the distribution of physical components of galax- 
ies in wavelength space. We find that it is valid to assume that the stellar and the gaseous 
components of galaxies span complementary subspaces. To first order, stellar and gaseous spec- 
tral features can be extracted from galaxy spectra through a simple matrix multiplication. By 
comparing the stellar continua obtained respectively using the model-based and the empirical 
approach through a commonality measure, we conclude that the latter may lose higher-order 
spectral features. 
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1. Introduction 

A common method of estimating physical parameters of galaxies is by least-square fit- 
ting stellar population synthesis models to their spectra. The approach is time-consuming 
when the model becomes large. This hinders the examination of all possible parameters 



defining the models (for example, PEGASE (Fioc & Rocca-Volmerange 1997) is defined 
by as much as ~10 parameters), especially for large sky surveys such as the SDSS which 
will produce 1 million galaxy spectra. On the other hand, principal components (PCs) 



are known to be powerful in classifying galaxy spectra (Connolly et al. 1995) This sug- 



gests that, as PCs form one kind of basis vectors, a subspace is a good description to 
the higher dimensional spectral space. In this work, we exploit the subspace methods in 
extracting physical information from galaxy spectra. 



2. Stellar and Gaseous Subspaces 

The projection matrix for a set of astronomical spectra (Yip et al. 2004) is defined 
to be P = J^e l e )( e l' where e are the eigenspectra (Connolly et al. 1995) of that set 
of spectra. We construct a projection matrix using the Bruzual & Chariot (2003) The 



subspace defined by this projection matrix, P(star), is referred as the stellar subspace. 
Each galaxy spectrum / in our sample (taken to be various types of SDSS galaxy spectra) 
is projected onto this stellar subspace, so that the resultant vector, /(star), is 



/(star) = P(star) x / 



(2.1) 



where x means a matrix multiplication. If the gaseous and stellar subspaces are comple- 
mentary, we have 

P(gas) = I — P(star) , (2.2) 
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Figure 1. The stellar (light gray) and gaseous 
(dark gray) components of the observed galaxy 
spectrum (black) obtained by the subspace 
method. The flux density of the stellar com- 
ponent is arbitrarily offset for clarity. 
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Figure 2. The commonality between the stel- 
lar subspaces obtained by different methods, 
as a function of the number of eigenspectra in 
each and all of the sets of spectra. 



where / is the identity matrix. The gaseous component can be obtained directly by 
/(gas) = F(gas) x /. 

Fig. Q] shows the resultant projected spectra /(star) and /(gas) for an example galaxy 
spectrum, which resemble the expected stellar and gaseous contributions. Since for a given 
subspace, its projection matrix is uniquely defined (Oja 1983) only a single projection 
matrix is needed to define an arbitrarily complex model. 



3. Definition of Stellar Subspace: Model-Based vs. Empirical 

As we construct the stellar subspace by adopting the stellar population synthesis mod- 
els, it is important to understand how well this model-based approach is when compared 
with other approach such as the empirical approach. Our chosen empirical approach is 
for a particular pixel by averaging the flux between the 40th and 60th quantiles of the 
flux distribution taken from a wavelength range of ±200 km/s. The emission lines are 
masked around Ao ±2000 km/s, with Ao being the wavelength of the line center. 

To measure the overlap between the stellar subspaces defined by different approaches, 
we use the commonality measure (see Yip et al. 2004 for its application) . Fig. shows 
that the least-square fit stellar continua (i.e. model-based) do not overlap 100% with the 
model spectra because of the sample variation arise from the selection of the galaxy. This 
will be our null measure, as such no other stellar continuum estimation method would 
exceed this line. The commonality between the empirical running-average approach for 
the stellar continua and the model drops below that of the null measure with the number 
of modes, indicating the loss of higher-order features in the former. 
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